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■ Abstract 



The results of tunnelling studies of the energy spectrum of two-dimensional 



> 

m 

, (2D) states in a surface quantum well in a semiconductor with inverted band 



\ structure are presented. The energy dependence of quasimomentum of the 

' 2D states over a wide energy range is obtained from the analysis of tunnelling 

^ \ conductivity oscillations in a quantizing magnetic field. The spin-orbit split- 

\ ting of the energy spectrum of 2D states, due to inversion asymmetry of the 

\ surface quantum well, and the broadening of 2D states at the energies, when 

- ^ , they are in resonance with the heavy hole valence band, are investigated in 

, structures with different strength of the surface quantum well. A quantitative 

analysis is carried out within the framework of the Kane model of the energy 
spectrum. The theoretical results are in good agreement with the tunnelling 
spectroscopy data. 

PACS number (s): 73.20. At, 73.40Gk, 73.40Qv 
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I. INTRODUCTION 



A specific feature of a semiconductor with inverted energy band structure is the absence 
of an energy gap between conduction and valence bands. These bands are the two branches 
of the fourfold Fg band. The states of Fg band are classified by the projection of the angular 
momentum J=3/2 onto the direction of quasimomentum k. The states with a projection 
of ±1/2 are conduction band states, hereafter referred to as the spin states of an electron. 
The Fg band, which is a conduction band in ordinary semiconductors A^B^ and A^B^ with 
an open gap, is the light hole band, and lies below the degeneracy point for the energy 
difference Eg = E^'' - E^». 

Such peculiarities of the band structure of inverted semiconductors lead to some special 
features of the energy spectrum of spatially confined systems, based on these semiconductors. 
Thus, it was predicted theoretically that so called interface 2D states can exist near the 
boundary of an inverted semiconductor even without an attractive electrostatic potential. It 
was shown in a number of theoretical and experimental articles that these states play 
a key part in forming the energy spectrum of 2D states in heterostructures and superlattices 
based on inverted semiconductors. 

The absence of a forbidden gap in the inverted semiconductors lead to the fact that at 
negative energies (hereafter we measure energy from the degeneracy point of Fg band in the 
volume) the 2D states in the surface quantum well are resonant with the heavy hole valence 
band and, therefore, broadening should be present in this energy range. Such effects have 
been discussed for narrow gap semiconductors. However in these materials, 2D electron 
states at energies less than the energy of the top of the valence band are separated from the 
heavy and light hole states by a barrier and the width of this barrier is proportional to the 
forbidden gap. So, the effect of the resonanse broadening in ordinary semiconductors should 
differ from that in gapless semiconductors. 

Of special interest are 2D systems in asymmetric quantum wells. The lack of inversion 
symmetry in this case leads to the spin splitting of the 2D subbands at A; 7^ due to spin-orbit 
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coupling even without a magnetic field. This phenomenon for 2D electron states has been 
studied extensively in semiconductors with Eg > 0. p|-p^ In such materials the conduction 



band is the twofold degenerate Fg band and the spin-orbit interaction can be taken into 
account within perturbation theory. |T^ As mentioned above, in inverted semiconductors 



the conduction band is one of branches of Fg band and the spin-orbit coupling has been 



considered from the onset. 15 



This work is devoted to the investigation of the energy spectrum of 2D states in the sur- 
face quantum well in the inverted semiconductor Hgi_^Cdi.Te using tunnelling spectroscopy 
in a quantizing magnetic field. This method was used for the first time in an investigation of 



the 2D states in InAs surface accumulation layer by Tsui in Ref. [16|. Since then, 2D states 



in a large variety of semiconductors, e.g. GaAs [|T^, InSb, |T8| and InAs, ||T9| were studied 



by tunnelling spectroscopy in a magnetic field. However there are only a few tunnelling 



experiments on structures based on inverted semiconductors. [|8|,p0 

Unlike traditional methods (galvanomagnetic phenomena, volt-capacitance experiments), 
which give information about the carriers at the Fermi energy, this method allows one to 
obtain information about the energy dependence of quasimomentum of carriers over a wide 
range of energies both for empty and occupied states. In this work we have also used a 
refinement of tunnelling spectroscopy which allowed us to investigate 2D states at energies, 
which are given by the applied bias, for different quantum well strengths. 



In the present article in contrast to previous publications PJ2^ we report on results 



obtained on tunnel structures with a larger strength of the surface potential well. The 
preparation of such tunnel structures has been made possible by the use of Yb as a metal 
electrode, which has a low work function. This allows us to observe both spin branches, 
which are split by the spin-orbit interaction, and 2D states at large negative energies where 
they are in resonance with the valence band. 

The article is organized as follows. In Section y some features of the tunnelling spec- 
troscopy of 2D states in a quantizing magnetic field, are considered. The experimental 
details and innovations in the tunnelling spectroscopy method are given in Section |ITT| . In 



Sec. experimental results obtained for the tunnel structures with different strengths of 
the surface potential well are presented and analysed. Section |V| is devoted to the theoretical 
description of the energy spectrum of 2D states in a surface quantum well of an inverted 
semiconductor. The basic equations, used in the calculations of resonant 2D states, are 
spelled out. A comparison between experiment and theory is discussed in Section |V1|, and, 



finally, conclusions are made in Section |VI1 



II. OSCILLATIONS OF THE TUNNELLING CONDUCTIVITY IN A 

MAGNETIC FIELD 

In principle, the bias dependence of tunnelling conductivity of a metal-insulator- 
semiconductor structure contains information on the energy spectrum of both 3D and 2D 
states which may exist in the surface quantum well of the semiconductor. The investigations 
of tunnelling conductivity oscillations in a magnetic field permit one to obtain the most re- 
liable information on the energy spectrum, because in this case it is possible to determine 
directly the positions of the Landau levels and energy intervals between them. 

Different types of oscillations of the tunnelling conductivity may occur in the metal- 
insulator-semiconductor structure with 2D states localized in the surface quantum well (Fig. 
|1]). The tunnelling current in such a structure is the sum of the current due to tunnelling 
to 2D (j2D) and to 3D (jsd) states. The oscillations of a^o = dj^n/dV versus magnetic 



field and bias have been considered earlier. [pi|-[23[| The maxima in a^n occur when the 



metal Fermi level is aligned with one of the bulk Landau levels of the semiconductor. These 
oscillations are periodic in the reciprocal magnetic field with the period 

where Ep is the semiconductor Fermi energy and S{Ep + eV) is the extreme cross-section of 
the iso-energy surface at the energy Bp + eV . In semiconductors with an isotropic spectrum 
S{E) = 7ck^{E). This gives one the opportunity to determine the energy dependence of the 



quasimomentum for the bulk states over a wide energy range. |2^j23[] 



The tunnelling current to 2D states at bias V, magnetic field B, and temperature T = 

is 

fEp+eV 

32d{V) = / g2D{E - Eo, B)D{E, V)dE (2) 

where g2D{E — Eq,B), Eq, D{E,V) stand for the density of 2D states, energy of the bottom 
of 2D subband and the barrier tunnelling transparency, respectively. An applied voltage not 
only shifts the Fermi quasilevels of the metal and semiconductor relative to each other but 
also changes the surface quantum well (and therefore Eq) and barrier transparency as well. 
Thus the tunnelling conductivity for tunnelling to 2D states is 

a2D = ^ = g2D{E - Eo, B)D{E, \/)U^+ey + 

+ Je. [ dV DiE,V)+g2oiE-Eo,B)^^jdE (3) 

To a first approximation, the variation of g2D{E — Eq, B) with bias is due to changes in Eq, 
and relationship between Eq and V is linear, i.e. Eo{V) = Eo{0) +aV and g2D{E — Eq, B) = 
g2D{E — EQiy), B). This results in the expression 

a2D = D{E, V)g2D{E - Eq{V), 5)^^+,^ + aD{V)g2DiE - Eq{V), 5)U^+ey - 

rEp+eV dD(E V) 

-aD{V)g2D{E ~ Eq{V),B)\ep + / g2D{E ~ Eq{V) , B) \ ; U e (4) 

JEp UV 

where D{V) is the mean value of D{E, V). A magnetic field B \\ n {n is the normal to the 
surface) quantizes the energy spectrum of 2D states and (?2D becomes an oscillatory function 
of B and E. One can see from expression (^) that two types of conductivity oscillation 
should occur in a magnetic field. 

(i) In the first type, oscillations arise whenever a Landau level of 2D states is aligned 
with the energy Ep + eV (the first and the second terms in Eq. (^). If the surface potential 
depends only slightly on V, then the positions of these oscillations allow one to immediately 
determine the energy spectrum of the 2D states by the same procedure as for 3D states (see 
Eq. ([I|)). As i? — s> the fan chart of these oscillations is extrapolated to the energy of the 
bottom of the 2D subband. 
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(ii) In the second type, oscillations appear whenever a Landau level of the 2D states is 
aligned with the energy Ep (the third term in Eq. (^)). As i? ^ the fan chart of these 
oscillations is extrapolated to the bias at which the 2D carriers (but not states) disappear, i.e. 
to the bias for which the bottom of the 2D subband becomes higher, than the semiconductor 
Fermi level. The amplitude of these oscillations is proportional to the rate of variation of Eq 
with bias (a). When a = 0, or i?o > Ep these oscillations are absent and only oscillations 
of the first type will be observed. 

We do not concern oscillations of the forth term in Eq. (^, because an oscillating function 
is an integrand, and, therefore, the amplitude of oscillations of this term is small. 

Up to now we have not considered possible oscillations of D{E^ V) in a magnetic field. 
The effects of these oscillations might be dominant in tunnel structures with a monocrystal 
barrier. For such a barrier the decay constant of the wave function in an insulator is rigidly 
determined by the energy. Therefore, small oscillations of the surface potential (ps (Fig. |l]) 
give rise to significant oscillations of D{E,V). Oscillations of as well as oscillations of 
the third term in Eq. (^ occur whenever a Landau level of the 2D states is aligned with the 
energy Ep. |T^ An analysis of the expression (^ shows that the second type of oscillations 
will be predominantly observed in such structures. We believe that the oscillations of this 
kind were observed in HgTe/HgCdTe heterostructures. [pi] 



III. EXPERIMENTAL DETAILS 

The differential conductivity ad = dj/dV as a function of bias and magnetic field in 
metal-insulator- Hgi_ xCdxTe (0.08 < x < 0.13) structures was investigated in magnetic 
fields up to 6 T at the temperature of 4.2 K. Tunnel junctions were fabricated on monocrys- 
talline p-Hgi_xCdxTe with the concentration of uncompensated acceptors Nj^ — Np, = 
(2.5 — 60) X 10^^ cm"^. The dopping level was determined from an analysis of galvano- 
magnetic phenomena in the temperature range of 1.5-70 K. Ultraviolet illumination for 5-15 
minutes in dry air was used to form a thin oxide, which served as a tunnelling barrier. Then 
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a metallic electrode (Yb) was evaporated through a mask. Several tunnel contacts prepared 
on each sample were investigated. The parameters of the tunnel structures are listed in the 
Table |. It is assumed that the surface electric field in these structures arises from the work 



function difference of Yb and the semiconductor. [25 



The resistance of our structures is 0.1-1 kOhm and is determined mainly by the bar- 
rier transparency because the resistance associated with transition of electrons between 2D 
and bulk states of the semiconductor is significantly less in metal-insulator-semiconductor 
structures based on inverted semiconductors. 

The traditional modulation procedure was used for measuring the differential conductiv- 
ity and its derivative da/dV . In some cases a refined method was used, In parallel with 
a small (about 1 mV) alternative voltage with frequency / = 670 Hz and direct bias V , im- 
pulses with a period T ^ 1// and duration t ^ 1// were imposed across the tunnel contact 
(Fig. These impulses lead to a change in the electric charge of the localized states which 
occur in the insulator or at the insulator-semiconductor boundary. If the relaxation time 
of these states far exceeds T, these impulses result in an increase or decrease (depending 
on sign of the impulses) of the strength of the surface potential well. Thus, measurements 
taken during the interval between impulses make it possible to get information about the 
energy spectrum of the 2D states with energy Ep + eV but for a different strength of the 
surface potential well. 



IV. EXPERIMENTAL RESULTS 

Let us consider first the results obtained from the measurements of the tunnel structures 
prepared on the most heavily doped sample (Table |). The oscillation curves for such 
structures are simpler and therefore easier to interpret. Typical magnetic field dependences 
of the oscillatory part of the derivative of the tunnelling conductivity with respect to voltage 
at different biases are presented in Fig. ^. Oscillations are observed for both orientations of 
magnetic field, B \\n and B n. The positions of the maxima of d'^j/dV'^ for i? || n in the 
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B, V coordinates are plotted in Fig. ^. It is clear that with B \\ n, every oscillation curve in 
Fig. ^ is a convolution of several components and therefore it is difficult to follow the maxima 
positions of each component in Fig. ^ For better distinction a Fourier transformation of the 
oscillation curves is carried out. As shown in Fig. |^ for i? ± there is only one fundamental 
field and for B \\ n two fundamental fields occur in the oscillation curves (a shoulder, which 
position coincides with the maximum position at -B _L n, is resolved at some biases). Now 
one can separate each component peak, take the inverse Fourier transform and follow the 
positions of the oscillation maxima of any one of the oscillation components independently. 
The curves in Fig. ^ were obtained in this way and it can be seen that they adequately 
describe all the experimental data. 

As mentioned above (Eq. (|l])) the period of oscillations A = Bj^ is determined by the 
value of quasimomentum. Bias dependences of k for both orientations of magnetic field are 
plotted in Fig. ^. 

For B J- n, the spectrum of 2D states is not quntized by the magnetic field and, therefore, 
the oscillations of the tunnelling conductivity are due to tunnelling to Landau levels of bulk 
states only. Hence open circles in Fig. ^ correspond to the energy spectrum of the bulk 
states. The experimental data are in good agreement with E{k) dependence calculated 
within the Kane model with parameters: P = 8 x lO^^eV cm, Aq = oo (where Aq is the 
energy difference between valence band Fg and split-off r7 band) and Eg listed in the Table 
|. Extrapolation of this curve to k = gives the semiconductor Fermi energy Ep = —eVo 
(Fig. P), which is also given in the Table |I[ 

For B \\ n, oscillations are connected mainly with tunnelling to 2D states localized in 
the surface quantum well of the semiconductor. This is evident from the typical angular 
dependence of the fundamental fields (Fig. |^). The reasons of some deviation from the 
classical angle dependence B{Q) = B{0)/cos{Q), where 9 is the angle between B and n, 
were discussed in Ref. [Q. At O > 20° the 6-independent maximum, which results from 
tunnelling to bulk Landau levels, is resolved in the Fourier transform (open circles in Fig. 

Thus, Fig. 1^ shows that two branches of 2D states exist in this structure and at A; — ^ 
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they are extrapolated to the bias Vq > 0. This means that there are no 2D states below the 
Fermi level, i.e. 2D electrons are absent and therefore only oscillations of the first type (see 
Section |T|) can be observed. 

There are two possibilities for the explanation of these two branches; they are the ground 
and excited subbands of 2D states or they are two branches of the ground subband split by 
the spin-orbit interaction in the asymmetric quantum well. An estimation of the difference 
between the energies of the ground and excited subbands at k = shows that for structures 
based on a heavily doped semiconductor with a small effective mass (see Table |) this value 



should be more than 50-100 meV (a comparison will be given in more detal in Section |V|). 
This value far exceeds the experimental one (Fig. Thus two branches of the 2D states 
in structure 10-1 (Fig. |^) correspond to two different spin branches of the ground subband 
which is split by spin-orbit interaction. Hereafter supercripts -|- and - refer to different spin 
branches. 

The structure 10-1 was also measured with voltage impulses which change the surface 
potential (see previous section) and shift the 2D branches relative to the bulk states. Fig. |^ 
shows this shift at V=80 mV versus the amplitude of the impulses. The distances between 
the branches of the 2D states and bulk states A{k^'^y = (k^'^)^ — k^uik versus quasimomen- 
tum of the lower branch of the 2D state (k^)"^ are plotted in Fig. The results for the 10-2 
to 10-6 structures, which differ by the value of the surface potential, are also shown in Fig. |^. 
It can be seen that when (A;"*")^ decreases as the result of reduction in the surface potential, 
the upper branch of the 2D states is pushed into the continuum, i.e. these localized 2D 
states disappear, and for the 10-3 to 10-6 structures, in which A{k~^y < 2 x 10^^ cm~^, the 
only one /c+-branch of 2D states exists. The bias dependence of the quasimomentum of 2D 



states for one of these structures (structure 10-4) over all bias range are shown in Fig. |T0. 

The more complicated oscillation curves for B \\ n were observed for tunnel structures 
prepared on lesser doped samples. The maxima positions in V, B coordinates are plotted in 
Fig. |TT]for structure 12-1 . An inspection of Fig. |TT] shows that two main oscillation types are 
observed. The maxima of these oscillations shift in opposite directions relative to V as the 
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magnetic field is varied. The behavior of the maxima positions with angle shows that both 
types of the oscillations are connected with 2D states. Taking the Fourier transform one can 
determine the fundamental fields, and using Eq. (ffl) one can calculate the quasimomentum 



of the states responsible for the oscillations at any bias (Fig. ^2|). The oscillations for 5 ± n 
as well as those for the structure 10-1 are due to tunnelling to the Landau levels of the 



bulk states. Thus the open circles in Fig. are the energy spectrum of the bulk states. 
The 1"*", 0~ and 0^ branches correspond to oscillations of the first type, i.e to oscillations 
arising whenever the Landau levels of the 2D states coincide with the energy Ep + eV . The 
branches a and h relate to oscillations of the second type, i.e. to oscillations which arise 
whenever Landau levels of 2D states coincide with the energy Ep, therefore branches O"*", 
h and 0~, a intersect at = (Fig. |l^). The decreasing of the quasimomentum at the 
energy E = Ep with V (branches a, h) is a result of the decreasing concentration of 2D 
electrons due to a reduction of the depth of the surface potential well as bias is increased 
(inset in Fig. |T2|). We reasone that two branches of the ground 2D subband, which is split by 
spin-orbit interaction, and one branch of the ecxited subband are observed for the structure 
12-1 (Fig. [T^). Such interpretation will be confirmed in Section |VJ The analogous results 
were obtained for the 13-7 structure. 

It is significant that at negative biases up to = —30 mV we observe oscillations of 
the first type which occur whenever the Landau levels of 2D states coincide with the energy 
Ep + eV. At these energies the 2D states in the surface quantum well of the inverted 
semiconductor are in resonance with the valence band states. However in spite of this 
fact broadening of the 2D states is not so large to destroy the oscillation picture. To our 
knowledge this is the first experimental observation of 2D states which lie significantly below 
the top of the valence band. 
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V. THEORETICAL MODEL 



To describe the energy spectrum of 2D system in wide-gap semiconductors, the one-band 
approximation is usually used. In this case the energy spectrum is parabolic. Such an ap- 
proach gives good results, because typical energies of 2D states in wide-gap systems are much 
less than the energy gap. In asymmetric quantum systems (e.g., in surface quantum wells), 
as a result of spin-orbit interaction, splitting of the energy spectrum arises even without an 
external magnetic field. It is common practice to interpret the experimental data in this 
case using the Bychkov-Rashba model. |^ The spin-orbit interaction is described here by one 
additional term, included in the dispersion law. It is linear in respect to quasimomentum, 
and contains a new parameter, which has to be calculated independently. 

The one-band model is inapplicable in the case of 2D states in narrow gap and inverted 
semiconductors. The strong interaction between the conduction and valence band makes it 
necessary to employ a multiband Hamiltonian in calculations of the energy spectrum of 2D 
systems. It is well known that the energy spectrum of A^B^ semiconductors can be described 
by the Kane Hamiltonian. Because the Dirac Hamiltonian, which is simpler for calculations, 
gives the same energy spectrum for electrons in the volume as the Kane Hamiltonian, the 
Dirac model is also widely used for calculations of the energy spectrum of 2D systems in such 
materials (see Refs. |T1|,|2^ and references therein). However, this is not quite correct. It was 
demonstrated in Ref. @] that parameters of the energy spectrum of 2D electrons obtained 
with the Dirac and Kane models, are radically different. The Kane model was used in our 
previous article on the energy spectrum of 2D electrons in HgTe/HgCdTe heterostructures. 



2^ There are a number of points to be made before applying this model in calculations of 



the energy spectrum of 2D system near the interface oxide/semiconductor. 

There is difficulty in choosing the boundary conditions at the oxide/semiconductor inter- 
face, because the energy band structure of oxide is unknown. The interaction with remote 
bands is usually neglected in the Kane model, when it is used for finding the spectrum of 
2D electron states. Then the Schrodinger equation is the system of ordinary differential 
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equations, and traditional boundary conditions (|\I/|=0 at the interface) lead in this approx- 
imation to the existence of unique solution: the wave function is equal to zero over all space. 

m 

One way out of this dilemma is to reduce the system of ordinary differential equations to 
one equation of second order which corresponds to one component of the wave function and 
to use zero boundary conditions only for this component. [|T0| , p9| -|3l| This is a good approx- 



imation for 2D electron states in wide gap semiconductors, in which one component of the 
wave function is much greater than the other. However, all the components are of the same 
order of magnitude in narrow gap and inverted semiconductors, and additional arguments 
are needed, as to why only one component has to go to zero at the semiconductor/oxide 
interface. 



Another approach has been suggested by Sobkowicz in Ref . p8| . Here, the assumption is 
made that the band structure of the insulator is similar to that of the semiconductor. The 
only difference is the value of the energy gap, which is much greater for an insulator, than 
that for a semiconductor. Moreover, the condition Dc, ^ E (where D^, are conduction 
and valence band offsets, respectively) seems to be natural, because neither electrons nor 



holes have to be emitted from semiconductor. Figure ITB| presents schematically the model 



band structure for the case of inverted semiconductors. Because the electrostatic potential 
is constant at z < and z > L, the exact wave functions are known in these regions. In this 
case the eigenvalue problem can be solved exactly using the techniques of direct numerical 



integration. [24 



Thus in the framework of this model one can understand, which insulator parameters 
correspond to the zero boundary condition for the second component of the wave function(see 
expressions (|^), which follow) used in Ref. ||31|. This can be understood from Fig. |T^, 



which shows the calculated dependence of the energy of the groung 2D subband on the 
value of Dj,, while the value of Dc is fixed (for details see Ref. [0]). The splitting of the 2D 



subband into two branches O"*" and results from spin-orbit interaction. It is clearly seen, 
that zero boundary condition used in Ref. |^ is a limiting case Di, oo and the limit is 
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reached very slowly. The results obtained in the above model come close to the solution of 
the zero boundary condition problem only at > 100 eV. 

The approach discussed above neglects the interaction with remote bands which corre- 
sponds to infinite heavy hole mass. This is suitable for calculation of the energy spectrum of 
2D electrons only for positive energies. When the energy of 2D states is negative, they are in 
resonance with the continuous spectrum of the heavy hole valence band. It is important to 
take remote bands into consideration in the calculation of the energy spectrum in this energy 
range, because the tunnelling of the carriers from the space charge layer to the volume of a 
semiconductor can lead both to a change in the energy and a broadening of the 2D states. 

To calculate the energy spectrum of 2D states taking remote bands into consideration, 
the usual assumption has been made, that the energy difference Aq between valence band 
Tg and split-off Fy band is very large. In this case the Kane Hamiltonian is a 6 x 6 matrix. 
We choose the direction z to be normal to the interface, and the direction of carrier motion 
along y. Then the Hamiltonian can be block diagonalized into two 3x3 Hamiltonians for 
two groups of states. The Hamiltonian matrix for the first group is defined by 



P (hi 9_ 

\2 dz 



\ 



dz) 



E^^ + eip{z)- 
-£71 {kl - 



92 

dp 







E^^ + eip{z)- 



V 



R2 



2m 



(5) 



where E^^, E^^ are the energies of corresponding band edges, 71 is parameter, which de- 
scribes the interaction with the remote bands, and (p{z) is electrostatic potential. The values 
P and 7i are assumed to be the same for both the semiconductor and insulator. The Hamil- 
tonian H~ for the second group of states is obtained from if + by replacing ky by —ky. Thus, 
the Schrodinger equation 
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(6) 



is a system of differential equations of second order, which determines two branches of the 
energy spectrum corresponding to two groups of states. The eigenvectors in the insulator 
and the volume of the semiconductor are known, because ip{z) = const at z < and z > L. 
There are two types of eigenvectors. The first type corresponds to a light particle and is 



where the energy for fixed /c^, and kz is obtained from the equation 

(e - E-^) [e - E^^ + ^-^{kl + kl)^ = \p\kl + kl). 

The second type corresponds to the heavy particle and is 

/ . \ 







VV3 



2p [k, + ''^^ 



for the energy 



(7) 



(8) 



(9) 



(10) 



The wave function in the insulator ^7 can be written as linear combination of ^^''^ and 



= ¥''\ky, -kf^) + C¥^\ky, -A;«) 



(11) 



where C is some as yet unknown multiplier. Here only the terms, which diminish at 2; — > 00, 
are given. 
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In the region z > L the wave function is given by 

^sc = C,¥''\ky, )) + C2^^'\ky, k^P) + 

+C3^^''\ky, + C^^^'^ky, -k^) (12) 



In Exps. (|rT]) and (|T^) k^ ^ k^J^^ stand for the quasimomentum components perpendicular to 
the interface, which for any given ky, and E are determined from (^) and (pi?D, respectively. 
Because only normalised solutions are of interest here, the coefficient C4 in (^) has to be 
zero. So, we can now state the problem which can be numerically calculated: at a fixed 
energy one needs to find a value of C for z < such, that subsequent numerical integration 
of the Schrodinger equation throughout the region of the potential results in a zero value 
for C4. 

The behaviour of the potential ^p{z) in semiconductor is determined from the Poisson 
equation for the charge density 

piz) = -e\NA - NolMiL - z) - 

-eE r gm'^iE,z)\'dE (13) 

J Ei 

where the summation runs over all occupied 2D subbands. Ei denotes the energy of the 
bottom of i-th subband. The second term describes the contribution of the electrons localized 
in the quantum well. In the absence of 2D electrons (this situation occures in the 10-1 to 
11-1 structures) the second term in ([l3|) is equal to zero, and the Poisson equation can be 
solved exactly. Then ^{z) is parabolic: 

\^s{l-z/L)\Q<z<L 
if{z) = I (14) 

I z> L 

where 

1 



eiNA-Nn), 
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and K is the dielectric constant. For the structures with a concentration of 2D electrons 
comparable to {N^ — N^) x L (as in structures 12-1 and 13-7) one needs to calculate the 



potential self-consistently. We use here the assumption made in Ref. Namely, at the 
calculations of the charge density distribution we suppose that the wave function is energy 
independent '^{E,z) = -^{Ep^z). 

Let us now consider peculiarities of the 2D states resulting from resonance with the heavy 
hole band. All the results demonstrated in this section has been obtained with the following 
parameter values Eg = —110 meV, P = 8 x 10~® eV cm, 71 = 2, A^^ — No = 1 x 10^*^ cm~^, 
K = 20, and a parabolic dependence of ^p{z) with ips = 275 mV. Dc = 2 eV and = 1 eV 
are used (the choice of these values will be justified in the next section). 

Figure |15| shows the z-dependence of |\l/(z)p for two energy values, corresponding to 
non-resonant {E > 0) and resonant {E < 0) 2D electron states. It is clearly shown that 
the wave function does not decay in the semiconductor region {z > L) at negative energies. 
This means that the charge carriers from the space charge layer may go into the volume 
of semiconductor, its wave function is transformed from an electron into a heavy hole wave 
function in the process. To find the broadening of the 2D energy levels associated with 
such a resonance one may turn to the scattering theory, namely Levinson's theorem. pB|,|52 



This allows a calculation of a density of states, added to the valence band by presence 
of the electrostatic potential ^p{z) and the semi conductor /insulator interface, treated as a 
scattering potential: 

I\'^{z)\'dz 

92D oc — (15) 

where A is the amplitude of at z > L and the integration runs over the space charge 
layer region. 



Figure |T6| shows the density of states corresponding to two groups of states, calculated 
for different values of quasimomentum k (in the following, k means ky). It can be clearly 
seen that the width of maxima rapidly increases with decreasing quasimomentum values. 
Peculiar feature is that the width of the maxima corresponding to the O"*" branch is greater, 
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than that of the branch at the same energy. This means that tunnelhng of 2D electrons 
to the valence band states is more effective for states of the 0^ branch. At k = the 2D 



energy levels become well defined again. 

The dispersion law E{k) of 2D electrons calculated for the described model with above 
parameters is presented in Fig. |T^. At E > the results coinside with those obtained without 
considering the remote bands. At negative energies when the 2D states are resonant, the 
energy levels are broad. This broadening is shown in Fig. as hatched regions. The width 
of these regions is defined here as the half-width of the peak of the density of 2D states (see 
Fig. |l^). In the narrow range of k"^ = (0.3 — 3) x 10^^ cm~^, the width of the density of 
states peak is very large, i.e. the 2D states are practically destroyed in this range. 



Earlier, the resonance broadening of 2D states was discussed in Ref. [28| for narrow 



gap semiconductors and in Ref. |^ for semiconductors with inverted spectrum. In these 
articles the tunnelling to light hole states was taken into account, but tunnelling to heave 
hole states was neglected. In such approximation the broadening of 2D states in inverted 
semiconductors occures only at the energy less than E^^ in the bulk of semiconductor. As 
shown above the tunnelling to the heavy hole states is effective and this process has to be 
taken into account in inverted semiconductors. 



VI. DISCUSSION 

Let us begin with the discussion of experimental results obtained for the tunnelling 
structures 10-1 to 11-1, based on heavily doped materials. 

As mentioned above there are two ways to understand the presence of two branches in 
the energy spectrum of the 2D states in the 10-1 structure. They are either ground and 
excited 2D subbands or two spin-split states of the ground 2D subband. The calculations 
carried out in the framework of the above model show (see Fig. |1^ that as the surface 
potential increases, the excited subband at A; = appears only when the energy of the 
ground subband is less than -100 meV (this value only shghtly depends on the values of 
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and Dy). This is in contradiction with the experimental results presented in Fig. ^. Indeed, 
extrapolating the experimental curves to A; = results in an energy difference between these 
brunches which does not exceed 10 meV. Thus, these two branches of the 2D states are two 
spin branches of the 2D subband split by spin-orbit interaction. 

For quantitative analysis of the experimental data the parameters D^, D^, and <fs{V) 
relationship should first be determined. The experimental results for the structures 10-1 
to 10-6 provide a way of estimating the parameters Dc and independent of (fs{V)- In 
reality, we can calculate the (fs relationships of the quasimomentum for the upper k~{(fs) 
and lower k~^{ips) branches of the 2D states at a fixed energy. Then eliminating ips we obtain 
a k~ versus dependence which is shown in the A{k~^'~y, {k~^y coordinates in Fig. |^ for 
different values of Df. and D^. The inspection of Fig. ^ shows that the results are not very 
sensitive to and as long as Df. and > Eg and ^ D^. But when > D^. and 
especially in the case when Dy = oo, which corresponds to the zero boundary condition for 
the second component of the wave function, the results of calculation significantly deviate 
from the experimental results. On the first sight it seems to be surprising that the authors 
of Ref. were able to explain the experimental results on the spin-orbit splitting of the 2D 
energy spectrum using such a type of the boundary condition. The analysis of this article 
shows that this is due to using of the approximate (quasiclassical) method of the solving 
of the eigenvalue problem in Ref. Indead, as is seen from Fig. ^ the results of the 

approximate calculations (dashed line) differ significantly from the exact solutions (line 4), 
but they are in a better agreement with the experimental data. 

Fig. ^ shows that the more suitable parameters are = 2 eV, Dy = 1 eV, and they 
will be used for the analysis of the results for all tunnel structures. 

We have calculated the bias dependence of the quasimomentum of the lower spin branch 
over the entire bias range for structure 10-1, using the surface potential as a fitting 
parameter (inset in Fig. With this 'PsiV) relationship we have calculated the position 
of the upper spin branch of the 2D states (Fig. One can see the good agreement with 
spin splitting over the whole energy range. Notice that at = the energy of the 2D states 
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coincides with the energy of the bottom of the conduction band, i.e. the binding energy of 
the 2D states at /c = is equal to zero. This is a specific feature of the 2D states in an 
inverted semiconductor for weak quantum well strengths. 

Curve 1 in Fig. |^ shows that the upper spin states must disappear, when A{k~^)'^ for 
the lower spin branch is less than 4 x 10^^ cm~^, which corresponds to a small surface 
potential (fs < 230 meV. Actually for the 10-3 to 10-6 structures only one spin branch exists 
over the whole energy range. This is illustrated by Fig. |10|, which shows the energy versus 
quasimomentum curve for one of these strictures (10-4). It can be clearly seen that the 
model employed describes the experimental data for these structures well. 

Let us inspect the results for the structure 12-1 (Fig. |12]). Contrary to structures of the 
first type (10-1 to 11-1) there are 2D electrons in this structure. Their concentration can be 
determined from the bias dependences of k"^ for oscillations of the second type (branches a 
and b in Fig. |1^. The concentration is varied from n2D = 1-3 x 10^^ cm~^ at V = —60 mV 
to = 1-0 X 10^^ cm~^ at = 50 mV. An estimation of the density of uncompensated 
acceptors in the space charge region gives a value of about 2 x 10^^ cm~^, which is of the 
same order of magnitude as n2D, therefore the self- consistency of the potential has to be 
included in the calculation. Knowing the density of 2D electrons at different biases (see 
section |^ and using a self-consistent procedure we have calculated the bias dependence of 
the surface potential (inset in Fig. |T2|). One can see that ips for the structure 12-1 is close 
to that of the structure 10-1 but due to the lower doping level (Table |) the width of the 
potential well is significantly larger so that the strength of the potential well is also larger 
and more than one 2D subbands are localized in such a well. Using this relationship between 
(fs and bias we have calculated the quasimomenta of 2D states over the entire energy range 



(curves 1^, 0~ and 0^ in Fig. |T^). Curves a and b in Fig. [T^ present the calculated values 
of the quasimomenta for both spin states at the energy E = Ep. As can be seen, there is 
good agreement for both types of oscillations over the whole bias range. Thus, both spin 
branches of the ground 2D subband and lower spin branch of the excited subband exist in 
the surface quantum well of this structure. 
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As was mentioned above (section 0), the 2D states should broaden at negative energy 



due to resonance with the valence band. As well as in Fig. |1^ this broadening is shown in 
Fig. |T2| by hatched regions. One can see that at the bias, where oscillations were observed 
experimentally, broadening does not exceed 5 meV, i.e. it is less than the cyclotron energy 
in the investigated semiconductors at the magnetic fields, where oscillations were detected 
{B > 1 T). With a more negative bias, the broadening becomes significantly larger and this 
is conceivable reason why we do not observe oscillations at these biases. 
Two comments are necessary in closing this section. 

At first sight it may seem that the model of the insulator used for calculations of the 
energy spectrum of 2D states in our structures is forced. What actually happens is that 
this model imposes some relationships between the components of the wave function at 
the interface. Any insulator considered in the framework of the kP-method imposes some 
relationships between the components too. The only point is then what relationships does 
it impose. In inverted and narrow gap semiconductors the amplitudes of all the components 
are about of the same value, therefore a widely used zero boundary condition for one of the 
component of the wave function |TI|pU| -pT[] seems to be artificial. 



All the experimental results on the energy spectrum of 2D states have been obtained from 
an analysis of oscillations of tunnelling conductivity in a magnetic field, but the calculations 
were carried out without a magnetic field. For the comparison of the calculated results 
with the experimental data, the assumption was used that the quasiclassical rule for the 
quantization of the energy spectrum of 2D states in a magnetic field was fulfilled. To check 
this assumption, calculations of the Landau levels were also carried out. It was shown that 
except for the low surface potential values, when the energy difference between 2D and bulk 
states with the same quasimomentum was small, the quasiclassical rule for quantization 
in magnetic fields was fulfilled with good accuracy. Some results of these calculations and 



analysis of the peculiarities of the behavior of 2D Landau levels were discussed in Ref. [20 
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VII. CONCLUSION 



The energy spectrum of 2D states localized in a surface quantum well in inverted semi- 
conductor was studied by tunnelling spectroscopy in a quantizing magnetic field. In the 
investigated structures the strength of potential well mainly depends on the width of the 
well which is controlled by the doping level: it increases when the acceptor concentration 
decreases. Two different cases were observed in our structures. 

In the structures based on heavily doped material only ground 2D subband exists and its 
bottom coincides with the bottom of the conduction band. Thus, in p-type material there 
are no 2D electrons and only empty 2D states exist. At 7^ the 2D ground subband is 
split by spin-orbit interaction. In structures with a larger potential well strength, both spin 
branches are observed, whereas in other structures the upper spin states are pushed into the 
continum and only the lower spin branch is observed. 

In the structures based on materials with A^^ — A^^i < 10^* cm~^ the strength of the 
potential well is so large that the ground and excited 2D subbands exist in these structures. 
In this case both spin branches of the ground subband and only lower spin branch of the 
excited subband are observed. In these structures the bottom of the ground subband lies 
significantly lower than the bottom of the conduction band and the 2D states are in resonance 
with the valence band over a wide energy range. The fact that oscillations relating to 
tunnelling to such states are observed, shows that the resonance broadening is not large and 
the 2D states are not destroyed at these energies. 

Theoretical calculations carried out in the framework of the Kane model describe well 
all the peculiarities of the energy spectrum and spin-orbit splitting of 2D states in inverted 
semiconductors. Theoretical calculations which take the finite value of the heavy hole effec- 
tive mass into account have shown that the broadening of 2D states due to resonance with 
the heavy hole valence band is large only over a narrow range of quasimomentum values. 
This is consistent with the results of our tunnelling experiments. 
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FIGURES 

FIG. 1. Energy diagram of a metal (M)-insulator-inverted semiconductor (SC) tunnel junction 
with a surface quantum well for a bias V > 0. 

FIG. 2. Voltage which is applied to a tunnel junction. The voltage is a sum of a dc bias V, an 
ac modulation voltage with amplitude Vm, and impulses with amplitude Vi. 

FIG. 3. d?j/dV'^ versus B curves at different biases for for the structure 10-1. 

FIG. 4. The positions of the maxima of the cPj/dV'^ versus B curves for the structure 10-1 for 
S II n. The curves are obtained in the way described in the text. 

FIG. 5. The results of Fourier transformation of the tunnelling curves plotted in Fig. ^. A 
shoulder on a low- field peak for a bias of 100 mV for 5 || n corresponds to tunnelling to the bulk 
Landau levels. 

FIG. 6. Bias dependences of A;^ for the structure 10-1. Open and full circles represent the 
experimental data for B J- n and B \\ n, respectively. The upper curve is the dispersion law of the 
bulk electrons, calculated in the framework of the Kane model with parameters Eg = —70 meV, 
P = 8 X 10^^ eV cm. The other curves are the result of theoretical calculations for the 2D states 
described in Section Inset shows the bias dependence of ifs- 

FIG. 7. Dependence of fundamental magnetic fields for the structure 10-1, measured at y = 80 
mV. Dashed line indicates the cos(0)~^ dependence. Open and full circles relate to tunnelling to 
the bulk and 2D states, respectively. 

FIG. 8. A quasimomentum of 2D states as a function of the amplitude of impulses for a bias 
of y = 80 mV, which corresponds to the energy E = Ep + eV = 65 meV, for the structure 10-1. 
Dashed line shows the value of k"^ of the bulk electrons for the same energy. Solid curves are merely 
guide for the eye. 
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FIG. 9. The differences A(A;"'")^ and A(/c~)^ as a function of (A;"*")^ for different structures, 
V = 80 mV. Points for structure 10-1 have been obtained at different amphtudes of impulses. 
Sohd curves are the result of numerical calculations carried out with parameters listed in the table 
and different band offset values: Dc = 2 eV, = 1 eV for curve 1, Dc = 5 eV, Dy = 1 eV for 
curve 2, Dc = 1 eV, = 1 eV for curve 3. The curve 4 and dashed curve are the results of exact 
and approximate (in accordance with Ref. |3^) calculations with zero boundary condition for the 
second component of the wave function. The upper curve is the same for the calculations with 
different parameters. 

FIG. 10. Bias dependences of k'^ for the structure 10-4. Open and full circles represent the 
experimental data for B J- n and B \\ n, respectively. The upper curve is the dispersion law of the 
bulk electrons, calculated in the framework of the Kane model with parameters Eg = —70 meV, 
P = 8 X 10~® eV cm. The other curve is the result of theoretical calculations for the 2D states 
described in Section Inset shows the bias dependence of ips- 

FIG. 11. Fan chart diagram for structure 12-1 at i? || n orientation. The solid curves are merely 
guide for the eye. 

FIG. 12. The bias dependences of k'^ for structure 12-1. Symbols show the experimental data: 
the open and full circles correspond to the bulk and 2D states respectively. The dotted curve 
is the dispersion law of the bulk electrons, calculated in the framework of the Kane model with 
parameters Eg = —110 meV, P = 8 x 10~^ eV cm. The solid curves are the theoretical bias 
dependences of k'^[Ep + eV) (curves 1"*", 0"*", 0^) and k'^{EF) (curves a, and h). Inset shows the 
bias dependence of (ps- The hatched regions show the broadening of the 2D energy levels. 

FIG. 13. Model of an insulator-inverted semiconductor structure with a surface quantum well 
used in the calculations. 
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FIG. 14. The dependence of the energy of two branches of the 2D ground subband at 
k = 3 X 10^ cm~^ on the valence band offset value while Dc is kept constant Dc = 1 eV. The 
dashed lines are the energies of the branches calculated under zero boundary condition for the 
second component of the wave function at the insulator /inverted semiconductor interface. The 
calculations have been carried out with parameters Eg = —70 meV, P = 8 x 10~^ eV cm, 
A^^ — Nd = 6 X 10^^ cm~^, ips = 240 meV, and for parabolic surface quantum well. 

FIG. 15. I^'p versus z curves for resonant (lower figure) and nonresonant (upper figure) 2D 
states. The parameters used are in the text. 

FIG. 16. The density of resonant 2D states for different quasimomentum values for 0"*" (a) and 
0~ (b) branches. 

FIG. 17. The dispersion law of the 2D states localized in a surface quantum well. Two subbands 

split by spin-orbit interaction are shown. The hatched regions show the broadening of the 2D energy 
levels. The 2D states of the 0"*" branch are distroied in the range 

/j2 = (0.3-3.0) X 10^2 cm due to 
high probability of tunnelling into the heavy hole states of the volume of inverted semiconductor. 

FIG. 18. The energy of the bottoms of the ground and first excited 2D subbands as a function 
of surface quantum well depth, calculated with parameters of structure 10-1. 
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TABLES 



TABLE I. Parameters of investigated structures. 



structure 


X 


Eg (meV) 


m„(102m)^ Na 


- Nd (lO^'^cm-3) 


Ef (meV) 




(IQ-Scm)'^ 


10-1 












275 ±5 


1.0 


10-2 












270 ±5 


1.0 


10-3 


0.125 


-70 ±5 


6.25 


60 


-15 ±5 


240 ±5 


0.95 


10-4 












230 ±5 


0.92 


10-5 












210 ±5 


0.88 


10-6 












180 ±5 


0.8 


11-1 


0.115 


-90 ±5 


8.0 


20 


-10 ±5 


200 ±5 


1.5 


12-1 


0.105 


-110 ±5 


9.75 


8 


-8 ±3 


275 ±5 


2.75 


13-7 


0.095 


-125 ±5 


11.5 


5 


-5 ±2 


260 ±5 


2.85 



'^The effective mass of electrons at the conduction band bottom in the volume of semiconductor, 
m is the free electron mass. 

''The values are given for V = 80 mV for the 10-1 to 11-1 structures, and for y = mV for other 
structures. 
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